_WEAPONS  ASSIGNMENTS  IN  A TIERED  AIMPOINT 
SYSTEM  WHICH  SATISFY  DISPERSION  CONSTRAINTS 
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RESEARCH  OBJECTIVES 

en  one  is  unable  to  place  the  available  weapons  to  cover  every 
installation  so  as  to  attain  a minimum  specified  damage  on  each  instal- 
lation, the  supply  of  weapons  may  be  considered  as  a scarce  resource 
which  must  be  assigned  to  targets  in  such  a way  as  to  optimize  a pre- 
scribed objective.  The  prescribed  objective  may  be  strictly  military  or 
economic  or  a composite  of  various  intentions.  On  the  other  hand,  when 
it  is  possible  to  sufficiently  fulfill  the  objective,  it  may  be  advan- 
tageous to  obtain  the  assignment  of  weapons  which  is  optimal  with  respect 
to  other  criteria  such  as  cost,  reserves,  and  combat  efficiency.  Within 
the  scope  of  this  study,  we  consider  a mixture  of  weapons  and  delivery 
vehicle  types. 

There  are  numerous  factors  and  constraints  which  pertain  to  this 
assignment  problem  in  its  actual  application,  and  such  actual  assignments 
require  quite  sizable  investments  of  manpower  and  computation  time.  This 
research  addresses  the  mathematical  modeling  analysis  and  solution  of  the 
assignment  problem  focusing  on  damage  evaluation  and  supply,  range,  and 
coverage  constraints  as  well  as  the  dispersion  restrictions  on  the  weapon 
carried  in  any  single  delivery  vehicle.  The  solution  of  the  model  is 
not  intended  as  a blueprint  for  actual  assignments  by  strategic  units,  but 
rather  by  concentrating  just  on  these  factors  we  intend  the  model  as  an 

analytical  tool  to  assess  the  present  detailed  assignment  algorithms 

. 

and  to  evaluate  future  policy  concerning  nuclear  weapon  and  delivery 
system  procurement,  deployment,  targeting,  and  strategic  arms  limitation. 
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Accordingly,  this  research  is  directed  toward  a model  which  affords  a 
reasonably  quick  solution  algorithm;  thus,  the  algorithm  may  be  used  to 
give  broad  consideration  to  many  hypothetical  alternatives,  to  provide 
a starting  point  from  which  a detailed  blueprint  for  actual  assignments 
by  strategic  units  is  constructed,  or  to  generally  evaluate  present 
assignment  schemes. 


STATUS  OF  RESEARCH 

At  the  beginning  of  this  research  the  problem  was  presented  in 
raw  form.  Although  many  important  components  had  been  identified,  these 
components  had  not  been  characterized  mathematically  and  there  was  no 
comprehensive  framework  for  synthesis  and  analysis.  Thus,  the  initial 
effort  was  concentrated  on  defining  and  mathematically  characterizing  the 
relevant  features  of  the  tiered  aimpoint  system,  weapons,  delivery 
vehicles,  and  a general  value  system.  These  characterizations  led  to  a 
preliminary  system  of  constraints. 

1.  Each  delivery  vehicle  carried  a fixed  number  of 
weapons . 

2.  The  weapons  could  only  be  assigned  to  DGZs  of  matching 
type  within  the  tiered  aimpoint  system. 

3.  Only  DGZs  in  range  of  the  delivery  vehicle  may  be 
assigned  weapons. 

4.  The  weapons  carried  by  any  single  delivery  vehicle  can 
only  be  disbursed  within  the  footprint  of  that  vehicle. 

The  consideration  of  range  resulted  in  the  Mercator  computation  of  acces- 
sibility of  weapon  to  DGZ.  While  the  footprint  of  any  delivery  vehicle 
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can  be  described  as  a complicated  amalgam  of  physical  processes,  such 
computations  have  proved  to  be  the  major  obstacle  in  providing  the 
desired  quick,  analytical  tool.  Herein  we  develop  the  parallelogram 
footprint  approximation  of  the  true  footprint;  this  device  allows  the 
consideration  of  footprint  in  the  model  without  excessive  computation- 
al burdens.  While  not  an  exact  description  of  reality,  the  parallelo- 
gram footprint  approximation  is  reasonable  for  the  purposes  of  the  quick, 
analytical  tool. 

Accordingly,  the  task  became  defined  as  finding  the  assignments  which 
obtain  the  optimal  or  near  optimal  level  of  value  subject  to  constraints. 
The  notion  of  value  extracted  due  to  the  assignment  of  any  particular 
weapon,  however,  was  obscured  by  the  possible  interactive  effects  of 
various  weapons . This  confounding  was  minimized  by  the  addition  of  dummy 
DGZs  to  resolve  the  cases  in  which  many  weapons  are  assigned  to  one  of  the 
original  DGZs  or  in  which  some  of  the  assigned  weapons  cover  the  same  in- 
stallations. While  this  adds  constraints  to  the  problem,  it  eliminates 
the  nonadditivity  of  the  values  extracted  by  each  weapon.  Dummy  DGZs 
allow  some  interactive  effects  among  weapons,  yet,  albeit  subtly,  allow 
the  separate  extracted  values  to  be  added  to  provide  a measure  of  utility 
of  an  entire  assignment. 

The  result  of  this  modeling  and  approximation  admitted  a mixed 
integer/linear  programming  mathematical  synopsis  of  the  problem.  While 

the  simpler  problem  without  dispersion  restrictions  can  be  solved  using 

branch-and-bound  with  linear  programming  as  the  fathoming  device  since 

hat  problem  structure  allows  quick  access  to  a feasible  integral  optimun, 

an  examination  of  some  small  sample  problems  revealed  that  this  algorithm 


would  be  slow  in  the  cases  including  dispersion  restrictions.  Problems 
of  this  type  contain  a large  number  of  variables.  Identification  of 
decompositions  of  the  problem  by  delivery  vehicles  enabled  the  construc- 
tion of  primal  and  dual  solution  algorithms  which  converge  to  the  optimal 
assignment.  So  even  if  the  solution  procedure  is  truncated  we  obtain  a 
feasible  assignment  and  a measure  of  its  closeness  to  optimality.  In 
view  of  the  analytical  objectives,  this  is  a reasonable  result. 


LIST  OF  WRITTEN  PUBLICATIONS  IN  TECHNICAL  JOURNALS 

"A  model  of  weapons  assignment  including  dispersion  constraints"  to  be 
submitted  to  Operations  Research  or  Naval  Research  Logistics 
Quarterly  is  in  preparation  by  Jerren  Gould. 
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INTRODUCTION 

There  are  various  types  of  delivery  vehicles  in  which  nuclear  weapons 
are  conveyed  to  a target.  Any  individual  delivery  vehicle  may  be  cate- 
gorized according  to: 

1 . launch  point 

2 . range 

3.  number  of  weapons  conveyed 

4.  yield  of  each  weapon 

5.  restrictions  on  the  dispersion  of  the  weapons 

Based  on  these  characteristics,  it  may  be  determined  whether  a particular 
collection  of  aimpoints  are  accessible  to  the  weapons  in  a given  delivery 
vehicle.  All  of  the  weapons  carried  in  a single  delivery  vehicle  must 
be  targeted  in  a neighborhood  of  the  others.  For  MIRVed  ballistic 
missiles  such  a neighborhood  is  called  a footprint;  for  manned  bombers 
a reasonable  flight  path  must  be  demonstrable.  For  delivery  vehicles 
which  carry  a single  weapon,  the  criterion  reduces  trivially  to  range. 

The  domain  of  a potential  adversary  contains  numerous  installations 
(airfields,  missile  silos,  factories,  etc.)  and  each  installation  has 
been  appointed  a value  according  to  a given  system  of  utility.  There  are 
several  levels  of  weapon  yield,  called  tiers.  The  value  extracted  from 
an  installation  depends  on  the  tier  of  the  weapon  assigned  to  cover  that 
installation  as  well  as  the  location  of  the  aimpoint  relative  to  the 
location  of  the  installation. 

The  main  treatment  herein  is  the  case  in  which  there  are  insufficient 
numbers  of  weapons  allocated  to  a strike  to  extract  value  from  every 
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installation.  Such  a scenario  is  called  a target-rich  environment.  The 
analysis  and  solution  of  other  important  scenarios  follow  analagously. 

TIERED  AIMPOINT  SYSTEM 

The  supply  of  available  weapons  consists  of  weapons  of  various  yields. 
The  tier  of  a weapon  corresponds  to  its  yield  inversely.  That  is,  the 
lowest  tier  corresponds  to  the  highest  yield  and  the  highest  tier  cor- 
responds to  the  lowest  tier;  we  shall  assume  that  the  tier  ranges  over 
' •»  • * the  integers  1 to  T . The  tiered  aimpoint  system  is  a device  to  cope  with 

the  diversity  of  weapon  yields.  For  each  tier  a collection  of  potential 
aimpoints,  or  designated  ground  zeros  (DGZ),  is  generated. 

The  domain  of  a potential  adversary  is  divided  into  complexes 
(target  islands) . A complex  is  a region  containing  installations  such 

t 

that  a hit  by  a weapon  at  any  point  within  the  region  will  not  affect  in- 
stallations in  any  other  complex.  The  location  of  tier  1 DGZs  are  deter- 
mined so  as  to  partition  the  collection  of  installations  in  the  complex 
into  sub complexes.  The  number  and  layout  of  these  DGZs  are  constructed 
so  as  to  obtain  at  least  a threshold  fraction  of  value  extracted  subject 
to  prudence  and  parsimony.  Within  each  sub complex  defined  by  the  tier  1 
DGZs,  the  higher  tier  DGZs  are  determined,  tier  by  tier,  in  the  same 
manner  as  the  tier  1 DGZs  are  found  within  the  complex.  The  value 
extracted  by  the  weapons  of  any  higher  tier  in  a subcomplex  must  approximate 
the  value  extracted  by  the  tier  1 weapon.  Various  parameters,  such  as 
hardness  of  any  installation,  affect  the  layout  of  the  DGZs  in  any  tier 
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and  heights  of  burst.  An  adaptation  of  Cooper's  algorithm  [13]  may  be 
used  to  define  the  layout  of  DGZs. 

Certain  hardened  installations  may  require  more  than  one  weapon  of 
the  same  tier  to  extract  sufficient  value.  In  such  a case  we  consider 
this  single  installation  as  several  installations.  Despite  the  fact  that 
these  dummy  installations  share  a common  location,  they  may  be  partitioned 
into  distinct  DGZs;  thus,  each  such  dummy  installation  may  be  covered  by 
a weapon  without  the  others  necessarily  being  covered.  The  values  ex- 
tracted from  these  dummy  installations  are  constructed  so  as  to  incorporate 
the  cooperative  effects  of  multiple  coverage.  For  example,  the  values 
extracted  from  the  additional  dummy  installations  may  represent  the  marginal 
value  extracted  due  to  the  assignment  of  the  additional  weapons  which 
cover  the  real  installation. 

We  label  all  of  the  DGZs  in  the  layout  as  j = 1,...,N.  Each  DGZ  has 
a unique  tier  even  though  the  coordinates  of  several  aimpoints  are  identical. 
Let  t2(j)  represent  the  tier  of  DGZ  j.  Let  I(j)  be  the  set  of  installations 
covered  by  DGZ  j.  By  construction,  for  each  tier  t,  (I(j)  | j € t2^(t) } is 
a partition  of  the  set  of  installations,  I.  Also,  if  t 2(3)  = 1 and 
j € I(j'),  then  I(j')  c I(j);  that  is,  the  higher  tier  DGZs  are  nested 
within  the  subcomplex. 

The  layout  of  DGZs  is  constructed  without  regard  to  the  delivery 
system  deployment , dispersion  characteristics,  and  inventory  of  weapons. 

If  there  is  no  assignment  of  weapons  to  DGZs  so  that  sufficient  value  is 
extracted  from  every  installation,  we  have  a target-rich  environment. 
Otherwise,  we  have  a weapon-rich  environment. 
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We  shall  assume  that  the  weapons  conveyed  within  each  delivery 

vehicle  are  identical.  We  label  the  delivery  vehicles  k = 1 K. 

Let  be  the  number  of  weapons,  or  reentry  vehicles  (RV) , carried  by 
delivery  vehicle  k,  and  let  t^(k)  represent  the  common  tier  of  weapons 
carried  by  delivery  vehicle  k. 

Each  DGZ  is  constructed  to  accept  a weapon  of  only  one  particular 
tier.  That  is,  a weapon  from  delivery  vehicle  k may  be  assigned  to 
DGZ  j only  if  t^(k)  = t2(j).  The  ranges  of  all  the  weapons  in  delivery 
vehicle  k are  identical.  Let  R^  be  the  common  range  of  weapon  in 
delivery  vehicle  k.  Since  the  locations  of  the  launch  site  of  delivery 
vehicle  k and  the  DGZ  j are  known,  by  spherical  trigonometry  and  considera- 
tion of  the  rotation  of  the  earth  during  flight  we  can  calculate  the 
distance  d^  between  the  launch  site  k and  DGZ  j.  It  is  a necessary  con- 
dition for  a weapon  to  be  assigned  from  delivery  vehicle  k to  DGZ  j that 
both  the  range  criterion  be  satisfied  and  the  tiers  of  the  weapon  and  DGZ 
match.  We  define  the  accessibility  index 

1 

“kj 


if  t^(k)  = t2(j)  and  d^^  _<  R^ 


otherwise. 

Clearly,  d^  need  only  be  calculated  for  delivery  vehicles  and  DGZs  with 
matching  tiers.  Another  computational  saving  may  be  accrued  whenever 
similar  delivery  vehicles  carrying  a common  tier  of  weapon  are  located  in 
a relatively  small  neighborhood  of  each  other  as  compared  to  the  flight 
distances;  in  that  case  we  can  say  that  the  delivery  vehicles  share  a 
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Figure  1 Mercator  Projection  with  the  North 
Pole  as  the  Launch  Site 


L 
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common  launch  site.  We  can  analogously  find  a representative  launch 


site  for  the  SLBMs  by  considering  the  centroid  of  any  patrol  region. 


If  one  wishes  to  ignore  time  of  flight,  the  computation  of 


may  be  simplified.  For  the  launch  site  of  each  delivery  vehicle,  or 


for  each  launch  field,  the  map  of  the  world  is  transformed  by  rectifying 


the  great  circles  through  the  launch  site  in  a manner  akin  to  a Mercator 


projection.  In  such  a case,  d^  is  readily  obtained  as  the  latitude 


coordinate  of  DGZ  j in  the  new  map  with  the  launch  site  as  a pole. 


VALUE 


If  a weapon  is  assigned  to  DGZ  j,  the  tier  of  the  weapon  must 


match  the  tier  of  the  DGZ.  By  construction,  vj^»  value  extracted 


from  installation  i when  a weapon  is  assigned  to  DGZ  j,  is  known. 


These  values  extracted  depend  on  the  value  of  the  installation,  its 


hardness,  and  its  geographic  relationship  to  the  aimpoint.  Random 


variations  due  to  aiming  and  reliability  uncertainties  are  incorporated 


by  fixing  v^  as  the  expected  value  extracted.  It  is  possible  that 


damage  may  be  incurred  to  installations  not  covered  by  the  DGZs,  but 


we  shall  assume  that,  in  terms  of  value  extracted,  this  damage  is  negli- 


gible; that  is,  v^  = 0 whenever  i i I(j).  Due  to  certain  strategic/ 


political  considerations,  certain  installations  must  not  be 
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: 


e 


damaged.  For  each  DGZ  j which  covers  such  a "flagged"  installation, 
we  fix  Vj^«0  to  insure  that  any  reasonable  weapons  assignment 
would  deny  weapons  to  any  DGZ  convering  installation  i. 

Let 


kj 


1 if  a weapon  from  delivery  vehicle  k is 
assigned  to  DGZ  j 

0 otherwise. 


Clearly,  the  decision  variables  a,  are  the  unknowns  we  must  determine 

Kj 

for  a solution  of  the  problem. 

The  values  extracted  from  distinct  installations  are  assumed  to 
be  additive.  However,  if  weapons  are  assigned  to  DGZs  which  cover 
the  same  installation  (because  of  dummy  DGZs  this  will  only  happen 
if  the  weapons  are  of  distinct  tiers),  the  total  value  extracted  from 
the  installation  should  not  be  the  sum  of  the  separate  extracted  values . 
Let  J(i)  = {j|i€I(j)J  be  the  set  of  DGZs  which  cover  the  installation 
i.  A reasonable  measure  of  the  value  extracted  from  installation  i is 
the  greatest  value  extracted  from  that  installation  for  any  weapon 
assigned  to  a DGZ  which  covers  installation  i;  that  is, 
max 


k.JS(i)  i1' 


Accordingly,  in  a target-rich  environment  we  wish  to  find  the  feasible 
assignment  {a^ } which  maximizes 


l max 
i k, j €J (i) 


"kJ  VJi- 


(1) 


Thus,  only  assignments  which  conduce  a largest  value  extracted  for  some 
Installation  contribute  to  the  objective;  hence,  others  should  be  reassigned 


*$•* ' 
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elsewhere.  Notice  that  the  a priori  creation  of  dummy  installations  for 
certain  targets  allows  some  value  extracted  to  accrue  when  the  target 
is  covered  more  than  once. 

Thus,  in  any  optimal  assignment  every  DGZ  which  receives  a weapon 
must  contain  an  installation  from  which  it  extracts  the  dominant  value. 
When  an  installation  is  hit  more  than  once,  according  to  the  objective 
(1)  some  weapons  effectively  extract  no  value  from  the  installation. 
Although  it  is  possible  for  this  case  to  occur  in  an  optimal  solution, 
in  a target-rich  environment  this  waste  should  rarely  occur.  That  is, 
almost  all  assignments  will  be  to  DGZs  in  which  there  is  no  competition 
among  the  weapons  to  extract  value  from  the  installations  covered. 

Hence,  as  an  approximation  we  exclude  such  competition  by  forbidding 
assignment  schemes  in  which  any  installation  is  covered  more  than 
once.  This  constraint  is 

I a < 1 Vi  (2) 

k,j€J(i) 

Accordingly,  the  objective  (1)  transforms  to 


I max  a v = £ a v , 

i k,j€J(i)  k,j  J 


(3) 


where  v.  = > v. . is  the  value  extracted  from  all  of  the  installations 

J i«(j)  J1 

in  DZG  j when  an  appropriately  tiered  weapon  is  so  assigned.  In  a 
target-rich  environment  we  may  presume  that  all  of  the  weapons  allocated 
to  a strike  will  be  used. 


14 


DISPERSION 


The  dispersion  problem  arises  because  the  weapons  from  a single 
delivery  vehicle  cannot  be  assigned  arbitrarily,  but  rather  each  must 
be  assigned  to  a DGZ  "near"  to  the  other  DGZs  receiving  a weapon  from 
the  same  delivery  vehicle. 

Each  weapon  carried  in  a MIRVed  ballistic  missile  is  first  con- 
veyed by  the  delivery  vehicle  until  it  is  near  the  target  region; 
each  weapon  separates  from  the  delivery  vehicle  and  is  finally  con- 
veyed to  the  assigned  DGZ  (aimpoint)  by  a reentry  vehicle.  The  region 
into  which  these  weapons  must  be  aimed  is  called  the  footprint  of 
the  delivery  vehicle. 

While  the  true  footprint  is  somewhat  amorphous  because  it  depends 
on  many  factors  which  cannot  and  should  not  be  considered  in  this 
analysis,  it  is  generally  somewhat  elliptical.  We  shall  assume  that  the 
boundary  of  the  footprint  is  truly  an  ellipse  with  an  axis  parallel 
to  a great  circle  through  the  launch  point.  The  approximation  of  the 
elliptical  footprint  by  a parallelogram  footprint  allows  a mathematically 
more  convenient  characterization  of  the  dispersion  constraints  for 
MIRVs.  The  parallelogram  footprint  allows  a linear  form  of  the  dis- 
persion constraints,  whereas  other  characterizations  of  the  footprint 
constraints  are  either  inherently  nonlinear  or  introduce  more  variables 
into  the  analysis. 


Parallelogram  footprint 

Suppose  delivery  vehicle  k 
of  its  footprint.  The  coordinates 


is  directed  to 

<VV 


(X^,Y^),  the  midpoint 


are  variables  to  be  determined 
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in  the  optimization.  Let  be  half  of  the  downrange  reach  and  p^/G^ 

be  half  of  the  cross  range  reach  of  the  parallelogram  footrpint  charac- 
terization of  delivery  vehicle  k.  Let  be  the  coordinates  of 

DGZ  j . DGZ  j will  be  within  the  parallelogram  footprint  if  and  only 


if  (X^)  is  chosen  so  that 

VvV  + 'v5*1  5 v 


Let  us  define 


lf  xr\ 

if  X <Xfc 


if  X ^ 


J tVxj  if  Xj<xk 

Clearly,  e k and  f are  nonnegative  and  not  simultaneously  positive. 
Clearly,  = lxj_Xk  I'  Similarly,  we  define 


y4-Yu  if  y >Y, 


« ,t<\ 


f 

tYk"yj 


if  yj5Yk 
if  *j<Yk 


and  observe  also  that  g and  h are  nonnegative,  not  simultaneously 

jk  j k 


positive,  and  W ■ 'vV- 
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Consider  the  following  constraints  to  represent  the  parallelogram 
footprint  constraints  for  the  delivery  vehicles 


*k  + ejk"fjk  “ xj 

Vj.k 

(7) 

Yk  + 8jk_hjk  = yk 

Vj,k 

(8) 

VW  + <W  - Vki  - V1"^  - 0 (9> 

where  is  a suitably  large  number.  It  can  be  shown  that  when  these 

constraints  are  embedded  within  the  total  model  that  there  are  algo- 
rithms which  can  guarantee  nonnegative  solutions  for  ejit»  ^jk’  ®jk’  ^jk 
with  the  appropriate  pairs  not  simultaneously  positive.  Hence,  the 
definitions  (5,6)  and  (7,8)  of  e , f , , g.,  , h , coincide.  Since 

JK  J«C  JK  JK 

a^j  shall  be  constrained  to  0 or  1,  these  constraints  (7-9)  admit 
two  cases: 

(1)  a = 1;  here  constraints  (7-9)  are  equivalent  to  the  fact 
that  DGZ  j must  lie  within  the  paralellogram  footprint  of 
delivery  vehicle  k as  it  must  if  such  an  assignment  is 

to  be  made. 

(2)  a.  = 0;  since  a weapon  from  delivery  vehicle  k is  not 

assigned  to  DGZ  j , it  is  irrelevant  whether  DGZ  j lies  in 
the  parallelogram  footprint.  Here  constraints  (7-9)  merely 
serve  to  reflect  the  fact  that  DGZ  j lies  within  a much 
larger  parallelogram.  So  if  is  chosen  sufficiently  large, 

these  constraints  will  always  be  satisfied  when  a^  **  0. 
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Finding  the  best  parallelogram  footprint  approximation  to  the 
true  footprint 

For  our  analysis  here  we  shall  consider  the  midpoint  0 as  the 
origin  of  a Cartesian  coordinate  system.  There  are  two  types  of  errors 
that  may  occur  in  the  designation  of  any  footprint  for  analytic  purposes, 
let  P denote  the  set  of  points  within  the  parallelogram  and  T 
denote  the  set  of  points  within  the  true  footprint.  We  have  a type  1 
error  for  the  points  in  T-P;  roughly  speaking,  a type  1 error  implies 
we  have  possibly  prohibited  making  an  assignment  because  of  the  mathe- 
matical formulation  when,  in  fact,  the  assignment  is  physically 
realizable.  On  the  other  hand,  we  have  a type  2 error  for  the  points 
in  P-T;  a type  2 error  implies  that  we  have  possibly  permitted  an 
assignment  which  is  not  physically  realizable. 


Figure  2 Comparison  of  true  footprint  and 
parallelogram  footprint  with 
coincident  midpoints 


k 


2 
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Let  be  the  loss  associated  with  a unit  area  of  type  1 error 

and  Cj  be  the  loss  associated  with  a unit  area  of  type  2 error.  Hence, 
a measure  of  closeness  of  the  parallelogram  approximation  is 

C(P)  = c^  x area  of  (T-P)  + Cj  x area  of  (P-T) . (10) 

Thus,  one  can  find  a best  parallelogram  footprint  approximation  by 

A A 

solving  for  P such  that  C(P)  £ C (P)  for  all  parallelograms  P with 
midpoint  0. 

Let  the  true  footprint  be  an  ellipse  with  half  downrange  reach 

r and  half  crossrange  reach  r /a.  Thus,  the  boundary  of  the  true 

2 2 2 2 

(elliptical)  footprint  is  a x +y  = r . The  boundary  of  a parallelo- 
gram footprint  with  half  downrange  reach  p and  half  crossrange  reach 

p/0  may  be  represented  in  the  first  quadrant  as  y=-0x+p.  Let 
22  22  22  1/2 

q = (ex  r +0  r -ct  r ) . The  points  of  intersection  in  Figure  3 are 

(x  .y^  = ' 21  2 (U) 

<7  +0* 


(x2*y2)  = ~2  2 

a +e 


(0p-q,CT  p+0q) 


t/»  ./• 


Figure  3 Quadrant  of  comparison  of  elliptic 
and  parallelogram  footprints 
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k. 


Hence,  C(p,0)  - 2c1t^“  (tan  *(”■)  “ tan  + ~q~  + X2P_  ^ ^ 


y,p 


y,p  2 , y , cx 

+ 2c2{x2p  — [tan-  (^)  + tan  {-+)]} 


(13) 


To  find  the  best  parallelogram  P(p,0)  we  must  solve  the  system 
of  equations 


3C(o.0) 

3& 


0 


dCjp^h.  = 0> 

30 


(14) 


Such  a solution  may  only  be  attempted  numerically  given  the  values  of 
Cl*  C2’  r*  an<*  a"  Accordlngly,  we  introduce  a simpler  alternative  method 
with  a different  measure  of  closeness. 

Here  we  restrict  the  problem  by  fixing  o=0;  that  is,  the  ratio 
of  major  axis  to  minor  axis  of  the  ellipse  and  the  parallelogram  are 
identical.  We  wish  to  find  the  value  of  p which  minimizes  the 
maximal  distance  disparity  between  the  boundaries  of  the  elliptical  and 
parallelogram  footprints.  By  disparity  we  mean  the  shortest  distance 
as  measured  from  any  fixed  point  on  one  footprint  boundary  to  the  other 
footprint  boundary.  By  symmetry  we  need  only  consider  the  first  quad- 
rant. From  Figure  4 it  becomes  clear  after  some  analysis  that  the 
maximum  is  obtained  as  either  the  length  of  AB,  CD,  or  EF.  From 
analytic  geometry  we  find 


length 

AB 

Vi  r-p 
(l+a2)1/2 

length 

CD 

p - r 

length 

EF 

(p-r)/o . 
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r/o  p/9 


Figure  4 Minimax  criterion  for  determining  the 
best  parallelogram  approximation  (o=9) 


Hence,  the  solution  for  p which  minimizes  the  maximum  distance  dis- 
parity reduces  to  the  following  linear  programming  problem: 

minimize  E (1 

subject  to  E > ^ \ wo  (1 

(14tj  ) 


E > p-r 
E > (p-r)/o 


r < p 


V2  r. 


The  value  of  p which  achieves  the  optimum  is 

~ = , [(14tj2)1/2  + V2  mind.a)  1 

p r 2 1/2 

[(1+ct  ) ' + min  (1,ct) ] 


So  according  to  this  criterion  the  best  parallelogram  footprint  approxi- 

2 2 2 2 

mat ion  to  the  true  footprint  a x +y  = r is  given  by  half  downrange 
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reach  p and  half  crossrange  reach  p/o.  In  applications  r and 
o depend  on  the  properties  of  delivery  vehicle  k;  hence,  p may 
depend  on  k.  In  considering  assignments  the  footprints  are  translated 
so  that  the  midpoint  of  the  footprint  of  delivery  vehicle  k is 

<W- 


MATHEMATICAL  STATEMENT  OF  THE  PROBLEM  FOR  BALLISTIC  MISSILES 


We  wish  to  find  assignments  of  weapons  to  DGZs  which  do  well 
(optimize)  the  total  value  extracted  subject  to  the  constraints  of 
accessibility,  avoidance  of  multiple  coverage  (with  the  a priori 
exceptions  due  to  dummy  DGZs),  and  dispersion  due  to  footprint.  Let 
denote  the  number  of  weapons  carried  by  delivery  vehicle  k; 
clearly,  we  must  also  include  the  constraint  (23)  that  no  more  than 
weapons  from  delivery  vehicle  k be  assigned  to  DGZs. 

Formally,  then  we  obtain  the  mixed  integer/linear  programming 
problem 


N K 

maximize  J £ a,  , v 
j=l  k=l  k;*  J 


subject  to  a^  < Vj,k 


(21) 

(22) 


N 


j=l  ^ 5 ^ 


Vk 


(23) 


K 

l l 


k=l  j€J(i) 


a <1 

akj  - 


Vi 


(24) 


X,  + en  ~f ” x. 
X jk  jk  j 


Vj.k 


(25) 
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\ + ®jk  * V ' 

Vj.k 

(26) 

8k<ejk  + £jk)  + <8jk  + hjk) 

-Vkj-  V1’*^  5 0 

Vj.k 

(27) 

Skj ,ejk*f jk,8jk’hjk 

Vj.k 

(28) 

a^  integer 

Vj.k 

(29) 

Since  is  constrained  to  be  a nonnegative  integer,  constraint  (24) 
further  implies  that  is  either  0 or  1,  as  is  desired.  It  becomes 
evident  now  that  there  exists  a feasible  assignment,  namely  a^  = 0 Vk,j; 
this  assignment,  however,  accrues  no  value.  The  objective  function  (21) 
drives  the  problem  to  find  the  maximal  total  value  extracted  subject  to 
the  constraints.  Since  the  set  of  feasible  assignments  is  nonempty  and 
finite,  problem  (21-29)  has  a solution. 

Due  to  the  large  number  of  variables  in  problem  (21-29)  and  the 
binary  nature  of  each  a^j,  the  computational  aspects  of  finding  the  solu- 
tion are  quite  cumbersome.  In  view  of  the  uses  of  the  results,  we  may  be 
satisfied  with  good  suboptimal  feasible  assignment  schemes.  This  approach 
is  especially  important  in  a target-rich  environment  in  which  a computer 
algorithm  may  spend  the  majority  of  time  finding  an  optimal  assignment 
with  small  marginal  return  over  a good  suboptimal  solution  reached 
quickly. 
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A DUAL  SOLUTION  METHOD 


Herein  we  shall  incorporate  the  accessibility  constraints  (22)  by 
removing  all  variables  subscripted  or  ^ whenever  = 0.  Let 

rk  = ^ ke  t*ie  nun*er  DGZs  accessible  to  the  weapons  of  delivery 

vehicle  k.  For  each  k we  may  define  the  functions  k(£) = (k,j^),  where 
is  the  £—  DGZ  accessible  to  delivery  vehicle  k according  to 
ordering.  We  now  introduce  the  following  notation: 


some 


4 = 

(ak(l) * * * ' » 

*L- 

(vk(i),  • . . , 

Vk(rk) 

(xk(l)’  * • ‘ » 

Xk(rk> 

4 = 

(yk(i)»  • • • » 

yk(rk) 

(Xk’  Yk’  ^(1)’ 

fk(l) ’ 

that 

the  components 

of  3c* 

’ ^(r^’  hk(rk)} 


component  of  k(£) . Suppose  that  m of  the  constraints  (24)  are  non- 


redundant.  Let  Qk  be  the  m x rk  matrix  with  element  (c,d)  being  1 if 


th 


DGZ  k(d)  appears  in  the  nonredundant  constraint  of  (24)  . Accordingly, 

the  mathematical  statement  may  be  rewritten  as: 
k 


■nax  I a. 
k=l 


subject  to  l Qk  ^ 

k=l 


(30) 

(31) 


J 
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'-H 


4 


For  A 


and  Z 


, the  constraints  (31,  32,  34) 


< j — i 

may  be  represented  as  FA  + HZ  = _b , where  ]>  = 


4/ 


and  A is 


binary.  Let  F = { A.|  3^2  0 3 HZ^  b - FA,  A binary}.  Hence,  A is  a 

feasible  binary  solution  iff  A € F. 

/x\ 


Let  n 


% 


Ac 


, where  is  m x 1 and 


4 

4 

>4 


with 


\ / 


a^,  .§k>  lxl,  r^  x 1,  and  2r^  x 1,  respectively,  for  k = 1,...,K. 

By  Farkas'  Lemma,  3Z^  2 0 ? HZ^  ~ b^  - FA  iff  (1>  - FA)  ' n 2 0 whenever 

H'n^  2 0,  2 0,  and  2 0,  2 0 for  k = 1, . . . , K.  Let  C be 


the  cone  of  such  values  of  n. 


By  duality  [1],  we  may  reformulate  the  problem  as 


K 


max 

ACF 


I v/  a,  + min  { (b^  — FA)  ' } 

k=l  n€C 


(35) 


C is  nonempty.  Moreover  since  a feasible  solution  to  the  original 
problem  exists,  the  minimum  in  (35)  is  obtained  at  an  extreme  point  of 


1 


C. 


Let  C have  extreme  points  ' , i « lf...,ne>  and  extreme  rays 

C r) 

, i = . Hence,  we  rewrite  (35)  as 


maximize  V 


subject  to 


V 2 J a£  + (b  - FA)'  4G>«  1 ° 1»-*-»ne 


(b  - FA)’  n^r)  > 0,  i - l,...,nr 
^ binary. 


Consider  now  a relaxation  of  (36-39) 


maximize  V 


subject  to 


V - l aj  v.  + (b  - FA)  » n.(e)  V i € I, 

k-1  ^ ‘ 


(b  - FA)  ’ n( ^ 2 0 V i € I2 
aR  binary, 

where  1^  and  I2  are  subsets  of  the  integers  l,...,ne  and  l,...,nr> 
respectively.  In  the  usual  way  of  considering  relaxations,  we  may 
obtain  an  iterative  procedure  for  solution  of  the  original  problem: 

1.  Initialize  1^  and  I2  with  a few  (or  no)  elements 
each. 

2.  Obtain  V°,  A^  as  the  finite  optimal  solution  to 
(40-43) . 
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Solve  W *»  min  (t>  - F^)  ' jn 
subject  to  H'ji  0 
\ 10,  4 > 0 Vk 

Ho  i 0 

(a)  If  W £ 0,  then  A^  is  the  solution  to  the  original 
problem  and  the  optimal  subjective  value  is  V°. 
Stop. 

(b)  If  0 < W < *,  then  (44-47)  has  a finite  optimal 


( c) 

solution  at  ji  ^ with 


(e) 


* X ^ **  + <-  • fv  ' * * 


B=1 

(e) 

The  vector  iv*  is  an  extreme  point  of  C not 

yet  considered  in  (41).  Enlarge  I,  so  that 

( c) 

n . is  considered  in  (41).  Return  to  Step  2. 


(44) 

(45) 

(46) 

(47) 


(48) 


(c)  If  W = «>,  then  (44-47)  has  an  extreme  ray 

(e) 

and  an  extreme  point  _n  such  that  the  objective 
becomes  unbounded  along  the  half-line 
n,  = + An(r),  A^O  . Accordingly, 

(b  - FA^) ' < 0.  Hence,  enlarge  I^  to 

include  . If,  in  addition,  (48)  is  true, 

(g) 

also  enlarge  1^  to  include  n.  . Return  to 
Step  2. 

The  finite  convergence  of  this  algorithm  follows  directly  from  the 
finite  number  of  constraints  in  (36-39). 
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We  may  simplify  problem  (44-47)  by  decomposing  it  into  sub- 
problems according  to  delivery  vehicle  since 

k=l 


\VV4A>\ 

k=l  k 


\V"Av  ■ '4 

k=l  k 


k*=l 


(49) 


is  separably  additive.  That  is,  solving  (44-47)  is  equivalent  to 
solving 

oy  - ' 4 + <4  ’!>  4 (50) 

k k 


subject  to 


(51) 

IV 

lo 

<* 

(52) 

min  (^“4’ak)ak 

(53) 

subject  to  > 0, 

(54) 

for  each  k=l K,  and 

K 

min  - ^Q^)’  Do 

(55) 

subject  to  r]^  > 0. 

The  quantity  W is  the  sum  of  the  optimal  values  of  the  objectives  of 
(50,  53,  55). 


A 
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Problem  (40-43)  is  a linear  program  with  binary  variables. 
Since  successive  problems  (40-43)  in  the  iteration  differ  by  the  addi- 
tion of  one  or  two  constraints,  a cutting-plane  method  [8,  11,  7,  9] 
may  be  efficient. 

At  any  iteration  V°  provides  an  upper  bound  on  the  objective 
value  of  the  original  problem.  As  we  proceed  through  the  iterations. 


V°  is  nonincreasing.  We,  however,  only  obtain  a feasible  solution  to 


the  original  problem  when  the  iteration  is  terminated. 


A PRIMAL  SOLUTION  METHOD 


Let  us  adopt  (30-34)  as  the  form  of  the  problem,  but  add 


the  redundant  constraints  a,  - J and  rewrite  (32)  as 

— k ^k 


, a,\  «.  /J 

B.  I “ T^kl  = c.  . 

k W = \h  1 -k 


Let  S, 


Iw 


have 


extreme  points 


fWa>' 


is  a bounded  convex  set;  an  arbitrary  element  can  be  written  as 


a convex  combination  of  its  extreme  points  ^£kY£k’  Wl,c,-c 

it 


where  Zy»,  =1 
l 


and  Y^-0.  Hence,  the  problem  without  the  integer  constraints  may 


be  stated  as 


V ■ maximize  Z Z Yg^  (P^ 


k l 


(57) 


30 


subject  to 


Z Zy-P,.  + I _ s 

. , k~  £k  mxm  — 

k l 


(58) 


zr 


Ik 


Y*k~° 


Vk 

v-e,  k, 


(59) 

(60) 


where  cp^  ■ i^k(a)  and  P^  = and  £ represents  slack 

variables.  The  Initial  basis  consists  of  s = J and  Y-,  i = * • 
corresponding  to  the  extreme  points  £^(a)  = — * Wit^ 

proper  choice  of  global  coordinate  system,  we  find  the  remaining  non- 
negative components  of  extreme  point  as  e ^ = x^  and  g^  = . 

The  form  of  the  master  program  (57-60)  is  readily  amenable  to  generalized 
lower  bounding,  a specialization  of  the  simplex  method  which  solves  such 
problems  while  maintaining  a mxm  "working  basis"  [5]. 

For  working  basis  8,  let  (w^  • u£)  = ^ B \ where 

consists  of  the  dual  variables  corresponding  to  the  m constraints  (58) 
and  consists  of  the  dual  variables  corresponding  to  the  K con- 

straints (59).  Let  j5^  denote  the  Kxl  vector  with  1 in  the  kttl 
row  and  0 otherwise.  Hence,  the  "added  value"  for  including  the 

variable  y.,  into  the  basis  is 
ok 


" ‘ 4 V L k(a)  • w. 


ok’ 


(61) 


Since  4.  .(a)  is  the  decision  variable  component  of  an  extreme  point 
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- 


of  S^,  we  may  find  the  largest  "added  value"  for  any  fixed  k by 


solving 


Wk  “ “x  K ~ -1  Qk^k  ' “ok 


subject  to 


a.  >0,  z,  ?0 
tc  — -^k  — 


(62) 


(63) 


(64) 


The  "added  value"  for  including  the  slack  variable  s^  into  the  basis 

is  -m  . If  max(max  W , max  (-w  ))<0,  then  we  have  obtained  the 
J k k j 3 

optimal  solution  to  the  problem  without  the  integer  constraints.  Other- 
wise, we  enter  into  the  basis  the  variable  corresponding  to  the  largest 
"added  value."  Note  that  the  solution  of  (62-64)  must  occur  at  an 
extreme  point  of  S^.  When  is  constrained  to  be  binary,  for  each  k 

it  can  be  arranged  that  there  is  a unique  in  constraint  (59) . 

For  each  k,  the  problems  (62-64)  are  much  smaller  than  the 
original  and,  hence,  the  solutions  are  computationally  easier  to  obtain. 
Therefore,  it  is  here  we  reintroduce  the  integer  constraints,  a^  binary. 
By  solving  (62-64)  with  the  additional  constraints  a^  binary  by  the 
method  of  cutting  planes,  we  adaptively  shrink  and  may  find  new 

"extreme  points"  due  to  the  cuts.  Note  that  it  is  not  necessary  in  the 
master  program  to  have  an  a priori  evaluation  of  all  the  extreme  points 
£ but  rather  to  keep  track  of  the  current  basis  and  entering 
extreme  point  as  defined  by  (62-64)  with  a^  binary. 

In  this  primal  algorithm  we  start  with  a basis  which  yields 
a 0 objective  value.  At  each  iteration,  the  objective  function  is 
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nondecreasing.  While  adding  the  constraints  a^  binary  at  each  iteration 
destroys  the  guarantee  of  convergence  to  the  optimum  of  the  original 
problem,  at  each  iteration  we  have  a feasible,  binary  solution. 


NEAR  OPTIMAL  SOLUTIONS 


At  any  iteration  in  the  dual  method  provides  an  upper 

bound  on  the  total  optimal  value  extracted;  however,  only  at  optimality 
is  the  corresponding  binary  A^  feasible.  At  any  iteration  in  the 
primal  method  V°  provides  a lower  bound  on  the  total  optimal  value 
extracted,  but  the  corresponding  A^  is  feasible  and  binary.  As 
the  dual  approaches  optimality  there  is  likely  to  be  a large  number 
of  iterations  with  marginal  reduction  of  total  value  extracted  due 
to  the  target-rich  environment.  Also,  in  addition  the  primal  method 
being  computationally  slower  because  of  the  combinatorially  greater 
binary  considerations  in  (62-64)  for  each  k,  there  is  no  guarantee 
of  convergence  of  the  primal  method  to  optimality.  Hence,  if  V° 
and  V°  are  relatively  close  we  may  terminate  computation  and  accept 
A^  as  a good  suboptimal  solution  with  total  value  extracted  V°. 
Alternatively,  it  may  be  possible  to  deassign  some  weapons  from  A^ 
to  obtain  a feasible  binary  solution  with  total  value  extracted  better 
than  V°  . The  degree  of  closeness  of  V°  and  V°  are  directly  related 
to  the  rates  of  convergence  and  the  computational  capability  and  patience 
of  the  user. 


ALTERNATIVE  PARALLELOGRAM  FOOTPRINTS 


In  the  model  (21-29)  previously  presented  we  have  restricted  each 
individual  delivery  vehicle  to  only  one  possible  footprint.  The 
physical  mechanism  whereby  weapons  are  thrust  from  delivery  vehicles, 
however,  allows  a variety  of  potential  footprint  shapes  ranging 
from  prolate  to  oblate.  Let 


m half  downrange  reach 
half  crossrange  reach' 


For  any  fixed  ct,  using  the  minimax  criterion  with  0 = ct  the  best 

parallelogram  footprint  approximation  may  be  represented  as  (0,p  ) 

0 


Figure  5 Oblate  and  Prolate 
Footprints 


MMM 
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As  we  shall  see,  each  possible  value  of  9 for  a delivery 
vehicle  introduces  a large  number  of  additional  variable.  Furthermore, 
these  footprints  include  much  of  the  same  area  in  common.  Let 
be  the  number  of  potential  footprints  for  delivery  vehicle  k.  From  the 
computational  standpoint,  the  problem  is  already  immense  with  *=  1 V k; 
however,  if  one  uses  vjc“2  for  some  k,  the  best  values  to  consider 
are  0^  (most  oblate)  and  9^  (m°st  prolate).  One  may  also  consider 
keeping  v ^=1 V k,  but  give  diversity  to  the  single  choices  of  6 ^ for 
identical  delivery  vehicle  types.  Recall  that  it  is  only  required  that 
the  relevant  DGZs  lie  within  the  parallelogram  footprint;  their  con- 
figuration within  the  footprint  is  irrelevant  for  the  purposes  of  this 


investigation. 

We  represent  the  m1"^  potential  parallelogram  footprint  as 

(Gj^.Prjjj)*  k=l,...,K  and  m=l,...,v^,  let  indicate 

whether  delivery  vehicle  k uses  its  m1*1  potential  parallelogram 

footprint.  Also,  let  d^m  indicate  whether  a weapon  from  delivery 

vehicle  k is  assigned  to  DGZj  using  footprint  m.  Note  that  a 

delivery  vehicle  may  adopt  only  one  of  the  potential  footprints.  Hence, 

we  may  generalize  the  model: 

K N 

maximize  Z Z a,  . v.  (6i 

k=l  j-1  J 


subject  to  a^  £ 


kjm  kj 


^kjm  “ ^km 
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N 

j-Aj  5 lv 


2 ^ a.  j - 1 

k=l  j€J(i)  KJ 


*k.  + eJkm  ' fllm  ' *1 


^km  + 8Jkm  ^jkm 


6l  (ej1  + ^11  ) + (g.l.  + !>,.  ) 

km  jkm  jkm  jkm  jkm 


-"kAj.  - s 0 


!,V  5 1 

m=l 


akj * ^kjm’  ejkm’  ^jkm’  ®jkm’  ^jkm’  ^km  ^ 


\j>  ’to-  ^kjm  lnte8er- 


1 


(69) 

(70) 

(71) 

(72) 

(73) 

(74) 

(75) 

(76) 

(77) 


FLIGHT  PATHS  FOR  MANNED  BOMBERS 

If  delivery  vehicle  k is  a manned  bomber,  then  it  proceeds 
from  a given  initial  point,  makes  a tour  of  its  assigned  DGZs,  and 
terminates  at  one  of  a set  of  designated  airfields.  If  we  allow  a 
manned  bomber  a generous  footprint  and  appropriately  calculate  a 
priori  accessibilities  we  may  defer  evaluation  of  feasibility  of  the 
flight  path  to  an  external  algorithm. 

< 
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If  the  terminus  Is  fixed  then  demonstration  of  flight  path 
feasibility  is  mathematically  equivalent  to  a traveling  salesman 
problem.  If  the  terminus  is  left  variable,  the  airfield  capacities 
must  be  observed.  While  such  algorithms  exist,  common  sense  is  likely 
to  be  quicker  because  ^ is  relatively  small  and  there  is  a natural 
order  to  the  tour.  Furthermore,  in  a target-rich  environment  we  may 
alter  the  assigned  DGZs  to  make  the  flight  path  feasible  or  more 
sensible  without  much  loss  of  total  value  extracted. 


1 
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NOTATION 


[ 


-k 


Vi 


4 

4 

c 

Vj 

dkjm 


jk 

ejkm 

F 

fjk 

^jkm 

8jk 

8jkm 

V 

h,. 
j km 

I 

I(j) 

I 

n 

J(i) 


indicator  of  the  DGZs  assigned  to  weapons  from  delivery 
vehicle  k 

indicator  whether  a weapon  from  delivery  vehicle  k is 
assigned  to  DGZ  j 

vector  of  assignment  variables  for  all  delivery  vehicles 
tentative  dual  assignment  solution 
tentative  primal  assignment  solution 

cone  of  dual  variables  which  yield  primal  feasible  solutions 

distance  between  the  launch  point  of  delivery  vehicle  k 
and  DGZ  j 


indicator  of  whether  a weapon  from  delivery  vehicle  k is 
assigned  to  DGZ  j using  footprint  m 

positive  part  of  x^  - X^ 

positive  part  of  x^  - X^m 

set  of  feasible  binary  assignments 

negative  part  of  x^  - X^ 

negative  part  of  x^  - X^ 

positive  part  of  y^  - 

positive  part  of  y - Y 

j km 

negative  part  of  y^  - Y^ 

negative  part  of  y - Y^ 

set  of  installations 

set  of  installations  covered  by  DGZ  j 

n x n identity  matrix 

set  of  DGZs  which  cover  installation  i 


J 


n x 1 vector  of  unities 


K 

* k(£) 

N 

n 

e 


P 


r 


T 


'ji 


number  of  delivery  vehicles 

a function  returning  k and  the  DGZ  accessible  to  weapon 
from  delivery  vehicle  k 

nunber  of  DGZs 

uumber  of  extreme  points  of  C 
number  of  extreme  rays  of  C 
region  enclosed  by  parallelogram  footprint 
coverage  induced  by 

installation  coverage  indicator  for  DGZs  accessible  to 
weapons  from  delivery  vehicle  k 

indicator  of  whether  delivery  vehicle  k adopts  potential 
parallelogram  footprint  m 

half  downrange  reach  of  true  footprint 

number  of  DGZs  accessible  to  weapons  from  delivery  vehicle  k 

range  of  weapon  from  delivery  vehicle  k 

polytope  generated  by  constraints  relative  to  delivery 
vehicle  k 

region  enclosed  by  true  footprint 

a function  such  that  t^(k)  is  the  tier  of  weapon  carried 
by  delivery  vehicle  k 

a function  such  that  t ^ ( j ) is  the  tier  of  DGZ  j 
a suitably  large  number 

value  extracted  from  the  installations  in  DGZ  j 

values  of  the  DGZs  accessible  to  weapons  from  delivery 
vehicle  k 


value  extracted  from  installation  i when  a weapon  is 
assigned  to  DGZ  j 


AO 


V°  tentative  dual  objective  value 

V°  tentative  primal  objective  value 


4 


crossrange  coordinates  of  DGZs  accessible  to  weapons  from 
delivery  vehicle  k 


<xrV 


coordinates  of  DGZ  j 


(Xj^.Y^)  midpoint  of  the  parallelogram  footprint  of  delivery  vehicle  k 


th 


()C  ,Y  ) midpoint  of  the  m^  potential  parallelogram  footprint  of 
delivery  vehicle  k 


4 

Z 


downrange  coordinates  of  DGZs  accessible  to  weapons  from 
delivery  vehicle  k 

vector  of  dispersion  variables  for  all  delivery  vehicles 


midpoint  of  the  parallelogram  footprint  of  delivery  vehicle  k 
and  differences  to  coordinates  of  DGZs  accessible  to  weapons 
from  delivery  vehicle  k 


dual  variable  corresponding  to  the  capacity  constraint  for 
delivery  vehicle  k 


“kj 


accessibility  indicator  for  weapons  from  delivery  vehicle  k 
to  DGZ  j 


4 


vector  of  dual  variables  corresponding  to  the  parallelogram 
footprint  constraints  for  delivery  vehicle  k 


vector  of  dual  variables  corresponding  to  constraints  of 
location  of  the  footprint  of  delivery  vehicle  k 


indicator  of  inclusion  of 
k^  column  of  1^ 
vector  of  dual  variables 

vector  of  dual  variables  corresponding  to  the  single  coverage 
constraints 


4 


vector  of  dual  variables  corresponding  to  capacity  and  dispersion 
constraints  for  delivery  vehicle  k 
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an  extreme  point  of  C 
C r) 

an  extreme  ray  of  C 

0^  ratio  of  dovmrange  reach  to  crossrange  reach  of  the  parallelo- 

gram footprint  delivery  vehicle  k 

0^  ratio  of  dovmrange  reach  to  crossrange  reach  of  m—  potential 

footprint  of  delivery  vehicle  k 

number  of  weapons  in  delivery  vehicle  k 

v.  number  of  potential  parallelogram  footprints  of  delivery 

vehicle  k 

311  extreme  point  of 

^£k^a)  elements  of  corresponding  to  the  assignment  variables  of 

delivery  vehicle  k 

j^k(z)  elements  of  corresponding  to  the  location  of  the  footprint 

of  delivery  vehicle  k 

p,  half  downrange  reach  of  the  parallelogram  footprint  of  delivery 

vehicle  k 

half  downrange  reach  of  the  m^  potential  parallelogram  footprint 
for  delivery  vehicle  k 

a ratio  of  downrange  reach  to  crossrange  reach  of  true  footprint 


value  accrued  at 

"added  value"  for  including 

dual  variables  corresponding  to  convexity  of 

dual  variables  corresponding  to  single  coverage  constraints 
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EPILOGUE 

The  goal  of  including  the  dispersion  constraints  while  simul- 
taneously finding  optimal  (or  near  optimal)  assignments  of  weapons  to 
DGZs  is,  indeed,  valuable.  This,  however,  introduces  a problem  of 
great  computational  size.  While  we  have  produced  here  a model  and 
solution  methodology,  it  is  yet  to  be  determined  whether  any  computer 
system  possesses  the  size  and  speed  to  generate  significant  results. 
Thus,  the  next  step  appears  to  be  the  coding  of  an  algorithm  for  execu- 
tion on  an  increasingly  large  sequence  of  test  problems.  In  such  a 
way  we  may  determine  the  magnitude  and  scope  of  problems  which  can 
be  approached  with  presently  available  computer  technology.  It  appears 
evident  that  every  computational  trick  and  device  must  be  explored  in 
order  to  place  the  larger  problems  on  a computer  and  have  them  executed 
in  reasonable  time;  efficient  coding  and  virtual  storage  are  identified 
as  key  elements  of  this  task. 

Even  if  the  problem  with  the  entire  set  of  weapons  and  DGZs 
allocated  to  a given  exercise  is  too  large  for  significant  computational 
results,  it  is  likely  that  a decomposition  to  smaller  subproblems  with 
a smart  recombination  will  produce  good  suboptimal  results.  Also,  if 
at  the  first  stage  we  identify  a good  assignment  without  dispersion, 
then  within  a target-rich  environment  it  may  not  be  difficult  to  perturb 
the  assignment  to  obtain  a good  assignment  which  also  satisfies  the 
dispersion  constraints.  That  is,  problem  (21-24)  is  amenable  to 
linear  programming  algorithms;  if  we  apply  branch-and-bound  using 


* 

3 


special  structure  of  the  polytope 
we  may  quickly  find  an  optimal  solution 
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